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We present precise measurements of CP and CPT symmetry based on the full dataset of ii" — > rnr decays 
collected by the KTeV experiment at FNAL. We measure the direct CP violation parameteri?e(e'/e) = (19.2 
±2.1) X 10-". We find the K^-Ks mass difi'erence Am = (5265 ± 10) X 10*5 fe-^ and the Ks lifetime = 
(89.62 ± 0.05) X 10"^^ s. We test CPT symmetry by finding the phase of the indirect CP violation parameter 
e, 4>t = (44.09 it 1.00)°, and the difi'erence of the relative phases between the CP violating and CP conserving 



decay amplitudes for — > tt^ 



/>+_) and for K it"it° (</>qo), A0 = (0.29 ± 0.31)°. These results are 



consistent with other experimental results and with CPT symmetry. 



1. Introduction 

Violation of CP symmetry occurs in the neutral 
kaon system in two different ways. The dominant ef- 
fect is the result of an asymmetry in the mixing of K'^ 
and if such that Kl and Ks are not CP eigenstates. 
This effect is parameterized by e and is called indi- 
rect CP violation. The other effect, called direct CP 
violation, occurs in the K tttt decay process and is 
parameterized by e'. Direct CP violation affects the 
decay rates of i^T ^ Tr+Tr" and K tt^tt^ differently, 
so it is possible to measure the level of direct CP vio- 
lation by comparing r]^ and t^oo- 
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Measurements of tttt phase shifts [Ij show that, in the 
absence of CPT violation, the phase of e' is approx- 
imately equal to that of e. Therefore, Re{e'/e) is a 
measure of direct CP violation and Im{e'/ 1) is a mea- 
sure of CPT violation. 

For small |e'/e|, Im{e'/e) is related to the phases of 
77+ _ and r/00 by 
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Experimental results have established that Re{e'/e) 
is non-zero 0, i, i, Q. In 2003, KTeV reported 
Re{e'/e) = (20.7±2.8) x 10"* based on data from 1996 
and 19970. We now report the final measurement of 
Re{e'/e) from KTeV. The measurement is based on 
85 million reconstructed K ^ tttt decays collected in 
1996 1997, and 1999. This full sample is two times 
larger than, and contains, the sample on which the 
previous results are based. We also present measure- 
ments of the kaon parameters Am and T5, and tests 
of CPT symmetry based on measurements of At/) and 
- (t>sw- 



For these results we have made significant improve- 
ments to the data analysis and the Monte Carlo simu- 
lation. The full dataset, including those data used in 
the previous analysis, has been reanalyzed using the 
improved reconstruction and simulation. These re- 
sults supersede the previously published results from 
KTeV 6]. In this presentation, we will focus primarily 
on improvements to the neutral mode analysis which 
have reduced the systematic uncertainty in Re{e'/e) 
relative to the previous KTeV result. 



2. The KTeV Experiment 

The measurement of Rele'/t) requires a source of 
Kl and Ks decays, and a detector to reconstruct 
the charged (7r"'"7r~) and neutral (tt^tt") final states. 
The strategy of the KTeV experiment is to pro- 
duce two identical K^ beams, and then to pass one 
of the beams through a "regenerator." The beam 
that passes through the regenerator is called the re- 
generator beam, and the other beam is called the 
vacuum beam. The regenerator creates a coherent 
\Kl) + P \Ks) state, where p, the regeneration ampli- 
tude, is a physical property of the regenerator. The 
regenerator is designed such that most of the if —> tttt 
decays downstream of the regenerator are from the 
Ks component. The charged spectrometer is the pri- 
mary detector for reconstructing K — > tt+tt" decays 
and the pure Cesium Iodide (Csl) calorimeter is used 
to reconstruct the four photons from K vr^Tr" de- 
cays. A Monte Carlo simulation is used to correct for 
the acceptance difference between K — > tttt decays in 
the two beams, which results from the very different 
Kl and Ks lifetimes. The measured quantities are the 
vacuum-to-regenerator "single ratios" for K tt^tt" 
and K tt^tt'^ decay rates. These single ratios are 

proportional to [77^ /pp and \r]Qo/p\'^, and the ratio 

of these two quantities gives i?e(e'/e) via Eq. [1] 



2.1 . The KTeV Detector 

The KTeV detector (Figured]) consists of a charged 
spectrometer to reconstruct K — > tt+tt^ decays, a 
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Figure 1: The KTeV Detector 



pure Csl electromagnetic calorimeter to reconstruct 
K 7r°7r° decays, a veto system to reduce back- 
ground, and a three-level trigger to select events. Two 
virtually identical neutral kaon beams are incident on 
the detector; a movable active regenerator is placed in 
one of these beams to provide a coherent mixture of 
Ki^ and Kg- In this manner, we collect Ki^ — > tttt 
and Ks tttt decays simultaneously so that many 
systematic effects cancel in the ratios used to calculate 
Re{e'/e). 

The KTeV spectrometer consists of four drift cham- 
bers and a large dipole magnet. It measures the mo- 
menta of charged particles with an average resolution 
of 0.4%. The K — tt+tt" reconstruction achieves a 
z-vertex resolution of 5-30 cm and a mass resolution 
of 1.5 MeV/c^. The Csl calorimeter measures the en- 
ergies and positions of photons from the electromag- 
netic decay of the neutral pions in K — > 7r°7r° decays. 
It has an average energy resolution of 0.6%. The re- 
constructed decay vertex of the neutral pion is directly 
related to the energies of the photons: 



Zj^n — Zcsl — 



ri2\/EiE2 



(3) 



The K ^ TT^Ti^ reconstruction achieves a z- vertex res- 
olution of 20-30 cm and a mass resolution of 1.5 
MeV/c2. 

2.2. Monte Carlo Simulation 

KTeV uses a Monte Carlo (MC) simulation to cal- 
culate the detector acceptance and to model back- 
ground to the signal modes. The very different 
and Ks lifetimes lead to very different z-vertex dis- 
tributions in the vacuum and regenerator beams. We 
determine the detector acceptance as a function of 
kaon decay vertex and energy including the effects 
of geometry, detector response, and resolutions. To 
help verify the accuracy of the MC simulation, we col- 
lect and study decay modes with approximately ten 



Figure 2: Data-MC comparison of fraction of energy in 
each of the 49 Csl crystals in an electron shower, (a) 
The fraction of energy in each of the 49 Csl crystals in 
an electron shower for data, (b) 2003 data/MC ratio (c) 
current data/MC ratio 



times higher statistics than the K — > tttt signal sam- 
ples, such as Kl — > ■K^e'^u and Kl ^ tt^tt'^tt^. 

Many improvements have been made to the MC 
simulation since the previous result was published in 
2003 [6]. We have improved the simulation of elec- 
tromagnetic showers to include the effects of incident 
particle angles and to simulate the effects of wrapping 
and shims in the Csl calorimeter. We have improved 
the tracing of charged particles through the detector 
with more complete treatments of ionization energy 
loss, Bremsstrahlung, delta rays, and hadronic inter- 
actions in the drift chambers. We have also have up- 
dated a number of parameters that go into the kaon 
propogation and decay calculations. 

The current Monte Carlo produces a significantly 
better simulation of electromagnetic showers in the 
Csl. Figure [2] shows the data-MC comparison of the 
fraction of energy in each of the 49 Csl crystals in 
a shower relative to the total reconstructed shower 
energy for electrons from Ki^ — > ir^e^v decays. The 
majority of the energy is deposited in the central crys- 
tal since the Moliere radius of Csl is 3.8 cm. These 
particular plots are made for 16-32 GeV electrons with 
incident angles of 20-30 mrad, but the quality of agree- 
ment is similar for other energies and angles. The 
data-MC disagreement improves from up to 15% for 
the 2003 MC to less than 5% for the current MC. 
This improvement in the modeling of electromagnetic 
shower shapes leads to important reductions in the 
systematic uncertainties associated with the recon- 
struction of photon showers from i^T — > tt^tt" decays. 
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Figure 3: tt+tt^ invariant mass distribution for K 
■K^-K~ candidate events. The data distribution is shown 
as dots, the K — > 7r"'"7r~(7) signal MC (MC Sig) is shown 
as dotted histogram and the sum of signal and background 
MC is shown as a solid histogram. 

3. Data Analysis 

The K — > TT+TT^ analysis consists primarily of the 
reconstruction of tracks in the spectrometer. The ver- 
tices and momenta of the tracks are used to calcu- 
late kinematic quantities describing the decay. The 
K — > TT+TT" invariant mass distributions for each 
beam are shown in Figure [31 

To reconstruct K ^ tt^tt^ decays, we measure the 
energies and positions of each cluster of energy in the 
CsL A number of corrections are then made to the 
measured particle energies based on our knowledge 
of the Csl performance and the reconstruction algo- 
rithm. The precision of the Csl energy and position 
reconstruction is crucial to the K tt^tt" analysis 
and has been improved significantly since the previ- 
ous publication. We use the cluster positions and en- 
ergies along with the known pion mass to determine 
which pair of photons is associated with which neu- 
tral pion from the kaon decay and calculate the decay 
vertex, the center of energy, and the tt^tt^ invariant 
mass. The K ^ tt^tt" invariant mass distributions for 
each beam are shown in Figure 31 

For K —^ tt^tt" decays, the z vertex is determined 
using only the positions and energies of the four pho- 
tons in the final state. Therefore, the measured z 
vertex is dependent upon the absolute energy scale 
of the Csl calorimeter. The energy scale is set using 
electrons from Kl — > ir^e^v decays. A small residual 
energy scale mismatch between data and Monte Carlo 
is removed by adjusting the energy scale in data such 
that the sharp edge in the z-vertex distribution at the 
regenerator matches between data and Monte Carlo 
as shown in Figure [5l The final energy scale adjust- 
ment for 1997 data is shown as a function of kaon 
energy in Figure [51 the average size of the correction 
is ^0.04%. As a result of improvements to the simu- 
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Figure 4: is' — > tt^tt" m^OjrO distributions for for data 
(dots) and signal MC (histogram) in the vacuum (left) and 
regenerator (right) beams. The sum of the background 
MC is also shown (stars). The dashed lines indicate our 
cuts. 
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Figure 5: Regenerator beam K z-vertex distribu- 

tion near the regenerator for 1999 data and Monte Carlo, 
(a) Uncorrected data, (b) Data with energy scale correc- 
tion applied. 



lation and reconstruction of clusters, the required en- 
ergy scale adjustment is smaller and less dependent on 
kaon energy for low kaon energies than in the previous 
analysis. 

Background to the K — > tttt signal modes is sim- 
ulated using the Monte Carlo, normalized to data 
outside the signal region, and subtracted. In this 
analysis, we use decays from coherently regener- 
ated kaons only; any kaons that scattter with non- 
zero angle in the regenerator are treated as back- 
ground. Scattering background is the same for both 
K ^ TT^7r~ and K tt^tt^ decays so it can be iden- 
tified using the reconstructed transverse momentum 
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Figure 6: Change in the final energy scale adjustment rel- 
ative to the 2003 analysis. The dashed line represents no 
data-MC mismatch. 



of the decay products in charged mode; we use 
K — > TT+TT^ decays to tune the simulation of scat- 
tering background on which we must rely in neu- 
tral mode. Non-TTTT background is present due to 
the misidentification of high branching-ratio decay 
modes such as — > TT^e^v, n^fx^i', and 



Kl vr'^vr'^Tr^. Background contributes less than 
0.1% ofK TT+TT- data and about 1% of K ^ tt^tt^ 
data. 

Table H] summarizes the systematic uncertainties on 
Re{e'/e). We describe the procedure for evaluating 
several important systematic uncertainties below. 



Source 


Error on 

K ^TT + 


Re{e'/e) (xlO^*) 


Trigger 


0.23 




0.20 


Csl cluster reconstruction 






0.75 


Track reconstruction 


0.22 






Selection efiiciency 


0.23 




0.34 


Apertures 


0.30 




0.48 


Acceptance 


0.57 




0.48 


Background 


0.20 




1.07 


MC statistics 


0.20 




0.25 


Total 


0.81 




1.55 


Fitting 


0.31 


Total 


1.78 



Table I Summary of systematic uncertainties in Re(e'/e). 

Acceptance: We use the Monte Carlo simulation to 
estimate the acceptance of the detector in momentum 
and z-vertex bins in each beam. We evaluate the qual- 
ity of this simulation by comparing energy-reweighted 
z-vertex distributions in the vacuum beam between 
data and Monte Carlo. We fit a line to the data- 
MC ratio of the z-vertex distributions and call the 
slope of this line, s, the acceptance "z-slope." A z- 
slope affects the value of i?e(e'/e) by producing a bias 
between the regenerator and vacuum beams because 
of the very different z-vertex distributions in the two 
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Figure 7: Comparison of the vacuum beam z distributions 
for data (dots) and MC (histogram). The data-to-MC 
ratios on the right are fit to a line, and the z-slopes (see 
text) are shown. All distributions are for the full data 
sample used in this analysis. 



beams; we use the known difference of the mean z val- 
ues for the vacuum and regenerator beams along with 
the measured z-slope to evaluate the systematic error 
on Re{e'/e). 

Figure [7| shows the measured z-slopes for the 
full Xl^tt+tt", Kl TT^e^v, Kl^tt^tt^, and 
Kl 7r°7r°7r° event samples. We use the tt+tt" z- 
slope to set the systematic uncertainty and measure 



the n^e^Ve z-slope as a crosscheck. For neutral mode, 
we use the high statistics tt'-'tt^tt^ mode to set the sys- 
tematic uncertainty because it has the same type of 
particles in the final state as tt^tt^ and is more sensitive 
than tt'^tt'^ to potential problems in the reconstruction 
due to close clusters, energy leakage at the Csl edges, 
and low photon energies. 

Energy Scale: The final energy scale adjustment 
ensures that the energy scale matches between data 
and MC at the regenerator edge, but we must check 
whether the data and MC energy scales remain 
matched for the full length of the decay volume. We 
check the energy scale at the downstream end of the 
decay region by studying the z-vertex distribution of 
tt'^tt^ pairs produced by hadronic interactions in the 
vacuum window in data and MC. To verify that this 
type of production has a comparable energy scale to 
K tt^tt", we also study the z- vertex distribution 
of hadronic tt^tt^ pairs produced in the regenerator. 
The data-MC comparisons of reconstructed z vertex 
for these samples are shown in Figure [S] 

To convert these shifts to an uncertainty in i?e(e'/e). 
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Figure 8: Energy scale tests at the regenerator and vac- 
uum window. The diilerence between the reconstructed z 
positions for data and MC is plotted for the K 
regenerator , and vacuum window tt^tt" samples. The 
solid point at the regenerator edge is the K tt^tt" sam- 
ple; there is no difference between data and MC by con- 
struction. The open point at the regenerator edge is the 
average shift of the regenerator tt^tt" samples for all three 
years. The points at the vacuum window are the shifts for 
the vacuum window samples for each year separately. The 
hatched region shows the range of data-MC shifts covered 
by the total systematic uncertainty from the energy scale. 
For reference, the data-MC shift at the vacuum window 
from the 2003 analysis is also plotted. 
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Figure 9: Comparisons of the reconstructed kaon mass vs 
2- vertex (top left), kaon energy (top right), minimum clus- 
ter separation (bottom left), and photon angle (bottom 
right) for 1999 data and MC. The values plotted are the 
difference between the reconstructed kaon mass for each 
bin and the nominal PDG kaon mass. 



we consider a linearly varying energy scale distortion 
such that no adjustment is made at the regenerator 
edge and the z shift at the vacuum window is that 
measured by the hadronic vacuum window sample. 
The average energy scale distortion we apply is shown 
by the hatched region in Figure [S] We rule out energy 
scale distortions that vary non-linearly as a function 
of z vertex as they introduce data-MC discrepancies in 
other distributions. The systematic error on Re{e'/e) 
due to uncertainties in the K tt^tt'^ energy scale is 
0.65 X 10^^; this is a factor of two smaller than in the 
previous analysis. 

Energy Non-linearity: Some reconstructed quanti- 
ties in the analysis do not depend on the Csl energy 
scale, but are sensitive to energy non-linearities. To 
evaluate the effect of energy non-linearities on the re- 
construction, we study the way the reconstructed kaon 
mass varies with reconstructed kaon energy, kaon z 
vertex, minimum cluster separation, and incident pho- 
ton angle. Data-MC comparisons for these distribu- 
tions for the 1999 data sample are shown in Figure 
[51 To measure any bias resulting from the nonlinear- 
ities that cause the small data-MC differences seen 
in these distributions, we investigate adjustments to 
the cluster energies that improve the agreement be- 
tween data and MC in the plot of reconstructed kaon 
mass vs kaon energy. We find that a 0.1%/100 GeV 
distortion produces the best data-MC agreement for 
the 1997 and 1999 datasets. Figure [TOl shows the im- 
provement in data-MC agreement with this distortion 
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Figure 10: Effect of 0.1%/100 GeV distortion on Mk vs 
E/f for 1999 data. The values plotted are the difference 
between the reconstructed kaon mass for each bin and the 
nominal PDG kaon mass. 



applied to 1999 data. The data-MC agreement in the 
reconstructed kaon mass as a function of kaon energy 
has been significantly improved compared to the pre- 
vious analysis in which a 0.7%/100 GeV distortion 
was required for 1997 data. 



4. Results 

The final KTeV measurement of i?e(e'/e) for the full 
1996, 1997, and 1999 combined dataset is: 

Re{e'/e) = [19.2 ± l.l{stat) ± l.8{syst)] x IQ-X'i) 
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Figure 11: Re{e'/e) in subsets of the data sample. Each 
point is statistically independent. The dashed line indi- 
cates the value of Re{t/t) for the full data sample. The 
97a run range includes the 1996 K tt^tt" data. 



i?e(e'/e), and Im[e'/e) in a single, z-binncd fit. The 
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Figure 13: Re{€'/e) in 10 GeV/c momentum bins. The 
dashed line indicates the value for the full data sample. 
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Figure 12: Re{e/e) consistency with beam intensity, re- 
generator position, magnet polarity, and track bend. The 
low and high intensity samples are from 1999 only and 
have average rates of ~ 1 x 10^^ protons/s and ~ 1.6 x 10^^ 
protons/s, respectively. Reg- left and reg-right refer to the 
position of the regenerator beam in the detector. Mag-l- 
and Mag- are the magnet polarity and in/out are the bend 
of the tracks in the magnet. In each of these subsets the 
K — > 7r°7r° sample is common to both fits; the errors are 
estimated by taking the quadrature difference with the er- 
ror for the full dataset. The dashed lines indicate the value 
of Re{e'/e) in the appropriate full data sample. 



= [19.2 ±2.1] X 10""*. 

(5) 

We perform several checks of our result by breaking 
the data into subsets and checking the consistency of 
the i?e(e'/e) resuh. To check for any time dependence, 
we break the data into 11 run ranges with roughly 
equal statistics. We divide the data in half based on 
beam intensity, regenerator position, magnet polarity, 
and direction in which the tracks bend in the magnet. 
We check for dependence of the result on kaon mo- 
mentum by breaking the data into twelve 10 GeV/c 
momentum bins. The Re{e'/e) results for these tests 
are shown in Figures [Til EH and [T31 We find consis- 
tent results in all of these subsamples. 

We also measure the kaon parameters T5, Am, (jj^, 



systematic uncertainties are evaluated using a proce- 
dure identical to that used for the Re{e'/e) measure- 
ment. CPT invariance is imposed a posteriori includ- 
ing the total errors of the parameters with their cor- 
relations to obtain a precise measurement of Am and 
T5. The results are: 

Am|cpt = [5269.9 ± 12.3] X 10-12 s, 

Ts Icpt = [89.623 ± 0.047] x lO^h/s, 

0+_ - [43.76 ±0.64]°, 

000 = [44.06 ±0.68]°, ^ ' 

d<p = 0, = [0.40 ±0.56]°, 

Acj) = -3/m(e7e) = [0.30 ±0.35]°. 
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